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During  many  clear  twilights,  much  of  the  solar  sky  is  dominated  by  pastel  purples.  This  purple  light’s 
red  component  has  long  been  ascribed  to  transmission  through  and  scattering  by  stratospheric  dust  and 
other  aerosols.  Clearly  the  vivid  purples  of  post-volcanic  twilights  are  related  to  increased  stratospheric 
aerosol  loading.  Yet  our  time-series  measurements  of  purple-light  spectra,  combined  with  radiative 
transfer  modeling  and  satellite  soundings,  indicate  that  background  stratospheric  aerosols  by  themselves 
do  not  redden  sunlight  enough  to  cause  the  purple  light’s  reds.  Furthermore,  scattering  and  extinction 
in  both  the  troposphere  and  the  stratosphere  are  needed  to  explain  most  purple  lights.  ©  2003  Optical 
Society  of  America 
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1.  Introduction 

Ever  since  the  Krakatoa  volcano’s  explosion  in  Au¬ 
gust  1883,  major  volcanic  eruptions  have  been  fol¬ 
lowed  by  reports  worldwide  of  extraordinarily  vivid, 
often  purplish,  skies  during  clear  twilights.  Here 
we  call  this  characteristic  feature  of  the  solar  sky 
during  post-eruption  twilights  the  volcanic  purple 
light.  Naturally,  volcanic  purple  lights  occurred 
long  before  the  Krakatoa  eruption,  and  scattered  ac¬ 
counts  of  these  date  from  at  least  the  early  16th 
century.  1 

After  the  Krakatoa  event,  19th-century  scientists 
quite  reasonably  speculated  that  the  eruption  in¬ 
jected  dust  into  the  upper  atmosphere  and  that  the 
twilight  Sun’s  light  was  scattered  by  the  resulting 
dust  layers  (or  cloud  droplets  with  dust  nuclei)  to 
form  the  purple  light. Some  early  authors  as¬ 
sumed  that  spectrally  selective  diffraction  by  these 
relatively  large  scatterers  caused  the  phenomenon’s 
red  component,®’^  whereas  others  viewed  the  scat¬ 
terers  primarily  as  a  means  for  scattering  already 
reddened  sunlight  toward  the  surface.®  Not  surpris¬ 
ingly,  still  other  writers  favored  a  combination  of  dif¬ 
fraction  and  single  scattering.® 
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Thus  began  the  long  and  close  association  of  the 
volcanic  purple  light  with  stratospheric  dust  and 
ash.'^-ii  Yet  even  during  volcanically  quiet  periods, 
pastel  purples  often  dominate  the  solar  sky  during 
clear  twilights,  and  stratospheric  dust  has  been  used 
to  explain  this  ordinary  purple  light  too.12-14  Al¬ 
though  some  authors  have  been  cautious  about  iden¬ 
tifying  these  stratospheric  scatterers  as  dust,®’i®  only 
fairly  recently  have  twilight  colors  been  linked  defin¬ 
itively  (but  not  exclusively)  to  scattering  by  sulfuric 
acid  droplets  in  the  stratosphere. 

2.  Phenomenology  of  the  Purple  Light 

Whatever  the  choice  of  scatterers,  all  existing  purple- 
light  models  emphasize  the  role  of  the  stratosphere, 
not  the  troposphere.  Even  authors  who  clearly  note 
the  troposphere’s  role  in  reddening  sunlight  tend  not 
to  dwell  on  it,  instead  understandably  focusing  their 
attention  on  the  geometry  of  scattering  in  the  strato¬ 
sphere. ®d4  One  of  the  few  exceptions  is  Deirmend- 
jian,  who  wrote  in  1973  about  the  “close  dependence 
of  twilight  features  on  the  optical  characteristics  of 
the  local  and  trans-horizon  troposphere,  and  their 
variation  not  only  along  the  vertical  but  along  the 
horizontal  direction  as  well,  and  .  .  .  the  highly  sub¬ 
jective  nature  of  individual  impressions  of  twilight 
colors,  relative  brightness,  and  their  changes  with 
time.”i® 

Indeed,  our  measurements  and  modeling  of  scores 
of  clear  twilights  indicate  that  both  the  troposphere 
and  stratosphere  are  crucial  to  explaining  the  forma¬ 
tion  and  the  absence  of  the  ordinary  purple  light. 
Furthermore,  the  subjectivity  that  Deirmendjian 
notes  necessarily  pervades  existing  qualitative  de¬ 
scriptions  of  the  purple  light.  This  subjectivity  has 
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made  purple-light  phenomenology  fairly  imprecise  to 
date,  and  this  imprecision  leads  to  confusion  about 
which  purple-light  features  are  common  and  which 
are  not.  A  corollary  problem  is  that  existing  descrip¬ 
tions  fail  to  mention  some  purple-light  phenomena 
that  we  observed  often  during  our  research. 

For  example,  Neuberger’s  diagram  and  description 
of  evening  purple  light  indicate  that  it  begins  as  a 
small  circle,  approximately  10°-15°  in  diameter 
above  the  Sun’s  azimuth.  Minnaert  makes  this 
spot’s  initial  view-elevation  angle  h  ~  25°  above  the 
astronomical  horizon^'*;  apparently  this  is  the  “purple 
spot”  variant  mentioned  by  Gruner  and  Kleinert.® 
However,  Minnaert  soon  dismisses  this  “purple 
patch”  as  the  least  common  of  five  ways  in  which  the 
purple  light  starts.  His  preferred  origins  include  an 
outgrowth  of  the  antitwilight  arch  that  “spreads  more 
or  less  invisibly  over  the  zenith  and  on  reaching  the 
west  becomes  visible  there”  (emphasis  added). 
Bafflement  is  an  understandable  reaction  to  such 
conflicting  (and  sometimes  invisible)  visual  descrip¬ 
tions.  In  fact,  only  once  during  our  research  did  we 
see  any  feature  that  could  be  described  as  a  “purple 
patch.”  It  occurred  well  after  the  onset  of  evening 
purple  light,  changed  shape  and  position  rapidly,  and 
clearly  was  caused  by  low-contrast  crepuscular  rays 
rotating  across  the  sky  (this  rotation  is  most  easily 
seen  in  time-lapse  photography). 

Adding  to  the  purple  light’s  phenomenological  con¬ 
fusion  is  the  older  practice  of  using  “intensity”  as  a 
synonym  for  “colorimetric  purity.”  Thus,  when  Neu- 
berger  and  Minnaert  refer  to  its  maximum  intensity 
occurring  at  unrefracted  solar  elevation  ~  —4°, 
they  actually  mean  that  its  purity  is  greatest 
then.  12.14  In  fact,  both  the  radiance  and  luminance 
of  skylight  decrease  nearly  logarithmically  during 
evening  twilight.  However,  quite  early  in  our  re¬ 
search  one  of  us  (Lee)  succumbed  to  this  long¬ 
standing  confusion  about  intensity  versus  purity  by 
replacing  “intensity”  with  “luminance”  rather  than 
“purity”  when  he  revised  Neuberger’s  purple-light  en¬ 
try  in  a  meteorological  dictionary.!® 

We  describe  the  evolution  and  shape  of  the  purple 
light  more  prosaically:  It  becomes  visible  impercep¬ 
tibly,  and  its  commonest  form  is  a  broad,  vaguely 
eyelid-shaped  region  centered  on  the  Sun’s  azimuth 
(although  attempting  to  assign  the  purple  light  any 
shape  is  problematic).  Figure  1  shows  the  typical 
color  range,  although  not  the  full  angular  extent,  of 
an  ordinary  purple  light.  Contrary  to  existing  de¬ 
scriptions  that  restrict  the  purple  light  to  an  azi¬ 
muthal  width  of  40°-80°  and  a  maximum  height  of 
—45°,  we  sometimes  observe  the  entire  twilight  sky  to 
be  purplish  (Fig.  2).  We  routinely  estimate  the  pur¬ 
ple  light’s  maximum  width  as  90°-150°,  and  we  can 
find  its  characteristic  colorimetric  signature  at  the 
zenith.  Equally  significant  in  our  research  is  the 
fact  that  the  purple  light  exhibits  great  diurnal  vari¬ 
ability,  sometimes  failing  to  appear  at  all  (Fig.  3). 
Time-lapse  photography  graphically  demonstrates 
one  feature  that  our  predecessors  commented  on,  the 


Fig.  1.  Photograph  of  an  ordinary  purple  light  taken  during 
evening  twilight  at  Chesapeake  Beach,  Maryland,  on  23  February 
1997  at  2250  UTC.  Unrefracted  solar  elevation,  ho  =  —0.39°; 
azimuth  relative  to  the  Sun,  4>rei  =  0°  fit  image  center. 


dramatic  setting  and  disappearance  of  the  purple 
light  as  evening  civil  twilight  ends. 

3.  Canonical  Purple-Light  Geometry,  Colorimetry,  and 
Photometry 

To  quantify  purple-light  colors,  from  1998  to  2001  we 
made  time-series  measurements  of  —100  clear  twi¬ 
lights’  spectra  at  four  sites:  (1)  Owings,  Maryland, 
(2)  the  U.S.  Naval  Academy  (USNA)  in  Annapolis, 
Maryland,  (3)  Spain’s  University  of  Granada,  and  (4) 
a  rural  site  near  Marion  Center,  Pennsylvania.!® 
Ah  sites  except  Owings  have  a  nearly  unobstructed 
2'7T-sr  field  of  view  (FOV)  of  the  clear  sky;  the  fisheye 
images  in  Figs.  2  and  3  show  the  USNA  site’s  view  of 
the  sky.  Table  1  provides  other  site  details.  Our 
spectroradiometers  were  two  Photo  Research  models 
PR-650  with  wavelength  ranges  of  380-780  nm  and  a 
spectral  step  size  of  4  nm.^® 

Because  the  ordinary  purple  light  can  be  found  at 
and  near  the  zenith,  we  sometimes  measured  hori¬ 
zontal  irradiances  with  a  cosine  receptor.  More  of¬ 
ten,  however,  we  measured  skylight  radiances  by 
aiming  the  PR-650  spectroradiometer’s  l°-FOV  tele¬ 
scope  at  the  solar  azimuth  (relative  azimuth,  c|)rei  = 
0°)  and  elevation  angle  h  =  20°.  This  choice  of  h  is 
arbitrary,  but  it  does  agree  with  existing  definitions 


446 


APPLIED  OPTICS  /  Vol.  42,  No.  3/20  January  2003 


Table  1.  Geographic  Details  of  Our  Measurement  Sites 


Site  Name 

Latitude 

Longitude 

Elevation  (m) 

Granada,  Spain 

37°  ll'N 

3°  35'W 

680 

Marion  Center,  Pa. 

40°  49'N 

79°  5'W 

451 

Owings,  Md. 

38°  41'N 

76°  35'W 

15 

USNA,  Md. 

38°  59'N 

76°  29'W 

18 

Fig.  2.  Fisheye  photograph  of  an  ordinary  purple  light  taken 
during  evening  twilight  at  the  USNA  in  Annapolis,  Maryland,  on 
6  October  1997  at  2257  UTC.  Solar  elevation,  /ig  =  —4.07°;  the 
photograph  is  centered  on  the  zenith. 


of  the  purple  light’s  location.  Using  Rozenberg’s  def¬ 
inition  of  twilight, 21  we  measured  skylight  spectra  at 
30-s  intervals  when  ^  -1-5°.  Depending  on  the 
value  of  /iq,  we  could  acquire  each  visible-wavelength 
spectrum  in  0.02-6  s. 


Fig.  3.  Fisheye  photograph  taken  during  evening  twilight  at  the 
USNA  on  2  October  1997  at  2302  UTC.  Solar  elevation,  /ig  = 
-3.85°;  the  photograph  is  centered  on  the  zenith.  In  contrast  to 
Fig.  2’s  twilight,  this  sky  has  no  discernible  purple  light.  Figures 
1-3  were  all  photographed  with  the  same  type  of  color  slide  film, 
Kodak  Elite  II  (ISO  100).  Exposure  times  in  Figs.  1-3  ranged 
from  0.008  to  0.5  s,  so  reciprocity  failure  and  its  possible  color  shifts 
are  unlikely  here. 


Figure  4  illustrates  the  purple  light’s  basic  scatter¬ 
ing  geometry.  In  it,  we  cut  through  the  Earth  to 
display  a  cross  section  of  an  atmosphere  with  greatly 
exaggerated  thickness  (in  reality,  the  stratosphere 
extends  from  altitudes  of  approximately  10-17  km  at 
the  tropopause  to  ~50  km  at  the  stratopause).  Par¬ 
allel  rays  of  sunlight  a-d  illuminate  a  radiometer’s  or 
our  own  line  of  sight  (LOS)  from  space  down  to  point 
t  above  the  surface  observing  site  s;  Fig.  4  does  not 
show  these  rays’  refraction.  Note  that  the  exagger¬ 
ated  atmospheric  thickness  in  Fig.  4  introduces  a 
minor  artifact:  Unlike  in  the  real  atmosphere,  none 
of  Fig.  4’s  tropospheric  rays  reenter  the  stratosphere 
before  they  reach  our  LOS.  Naturally,  radiances  re¬ 
ceived  at  the  surface  are  the  result  of  extinction  losses 
and  scattering  gains  all  along  the  given  LOS. 

In  Fig.  5  we  plot  two  typical  chromaticity  curves 
that  result  from  our  time-series  spectra  of  ordinary, 
or  non-volcanic,  evening  twilights.  Although  we  use 
the  perceptually  isotropic  CIE  1976  uniform- 
chromaticity-scale  (UCS)  diagram,  ordinate  and  ab¬ 
scissa  scales  differ  in  Eig.  5  and  subsequent 
chromaticity  diagrams  in  order  to  show  as  much  de¬ 
tail  as  possible. 22  Line  markers  such  as  the  X’s  or 
small  circles  in  Eig.  5  and  other  figures  distinguish 
the  chromaticity  curves  from  one  another  but,  unless 
indicated,  have  no  other  significance.  Some  useful 
colorimetric  landmarks  in  Eig.  5  are  the  curved 


radiometer’s 


Ray  a  is  tangent  to  earth  at  this  hg, 
ray  c  is  tangent  to  tropopause. 


stratosphere 


troposphere 


s 


radiometer 


Fig.  4.  Simplified  scattering  geometry  for  the  purple  light;  atmo¬ 
spheric  thickness  and  Sun  elevation  are  exaggerated  here.  Before 
they  reach  the  radiometer’s  line  of  sight,  parallel  rays  of  sunlight 
are  reddened  to  var}ring  degrees  by  transmission  through  the  tro¬ 
posphere  or  stratosphere. 
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Fig.  5.  Portion  of  the  CIE  1976  UCS  diagram,  showing  temporal 
trends  in  purple-light  chromaticities  measured  during  evening  twi¬ 
lights  at  Owings,  Maryland  (29  December  1999),  and  Granada, 
Spain  (7  June  2000).  Each  spectroradiometer’s  FOV  is  1°,  and 
throughout  these  measurements  its  view  elevation  is  h  =  20°  and 
its  relative  azimuth  is  4),gi  =  0°.  The  cross  labeled  “ET  sunlight” 
is  the  chromaticity  of  sunlight  measured  above  the  atmosphere. 


Planckian  locus,  two  color  temperatures,  the  color  of 
extraterrestrial  sunlight  above  the  atmosphere  (de¬ 
noted  “ET  sunlight”),  and  a  horizontal  line  labeled 
“1  JND.”  This  last  is  a  perceptual  ruler  whose 
length  is  the  mean  MacAdam  just-noticeable  differ¬ 
ence  (JND)  along  our  chromaticity  curves. 
Throughout  this  paper,  the  scales  of  JND  lines  are 
those  of  the  accompanying  u'  axis. 

The  fairly  vivid  Owings  purple  light  of  12-29-99 
begins  at  local  sunset  (ho  =  -0.75°)  as  a  sharp  re¬ 
versal  of  the  trend  that  makes  skylight  increasingly 
blue  for  Hq  <  +1.18°.  These  purple-light  chroma¬ 
ticities  cross  the  Planckian  locus  en  route  to  maxi¬ 
mum  vividness  (i.e.,  maximum  colorimetric  purity)  at 
ho  =  —4.37°  and  then  recross  it  soon  after  civil  twi¬ 
light’s  end  {ho - 6°).  In  fact,  this  pattern  typifies 

the  evolution  of  a  vivid  ordinary  purple  light:  (1)  as 
the  Sun  nears  the  astronomical  horizon,  skylight  tens 
of  degrees  above  the  horizon  reddens  until  ho  ~  1°- 
2°;  then  (2)  this  skylight  gets  bluer  until  sunset  or  a 
few  minutes  thereafter  (<5  min),  followed  immedi¬ 
ately  by  (3)  the  start  of  a  purple  light  whose  (4)  purity 
maximum  occurs  at  ho  ~  -3.9°  (see  Table  1)  and 
ended  by  (5)  skylight’s  rapid  bluing  through  the  end 
of  evening  civil  twilight. 

Now  we  must  reconsider  the  term  “purity.”  In  our 
observations  the  apparent  purity  maximum  for  any 
given  purple  light  occurs  when  its  post-sunset  chro¬ 
maticities  reach  their  maximum  CIE  1976  u'  (or  1931 
x)  value.  In  other  words,  the  best  purples  are  seen 
at  the  post-sunset  “nose,”  or  rightward-most  excur¬ 
sion  of  the  chromaticity  curve.  We  define  this  ex¬ 
treme  point  on  a  twilight  chromaticity  curve  as  its 
purple-light  maximum.  During  the  12-29-99  Ow- 


Fig.  6.  Purple-light  luminance  and  integrated  radiance  as  func¬ 
tions  of  ho  for  the  29  December  1999  Owings  evening  twilight 
whose  chromaticities  are  shown  in  Fig.  5. 


ings  twilight,  this  occurred  at  ho  =  -4.37°,  and  it  is 
marked  with  an  arrow  in  Fig.  5.  Even  though  the 
corresponding  u' ,  v'  chromaticity  is  closer  than  oth¬ 
ers  to  such  nominal  achromatic  points  as  Fig.  5’s  ET 
sunlight,  the  maximum  u'  chromaticity  is  consis¬ 
tently  the  purple  light’s  most  visually  striking  color. 
One  can  easily  (and  legitimately)  reconcile  observa¬ 
tion  and  colorimetric  terminology  simply  by  picking 
another  achromatic  point  (e.g.,  use  u'  =  0.19  and  v'  = 
0.44  in  Fig.  5).  The  chromaticity  of  this  perceived 
achromatic  stimulus  itself  depends  on  surround  col¬ 
ors  or  on  those  of  the  ambient  illumination  (i.e.,  hor¬ 
izontal  irradiances).  Given  the  large,  slowly 
evolving  color  changes  that  occur  in  twilight’s  ambi¬ 
ent  illumination,  local  chromaticity  extrema  such  as 
Fig.  5’s  purple-light  maxima  are  easy  to  interpret  as 
fairly  saturated  colors,  whatever  their  absolute  chro¬ 
maticities.  This  apparent  conundrum  about  colori¬ 
metric  purity  demonstrates  that,  within  reasonable 
limits,  the  temporal  chromaticity  curve’s  shape  mat¬ 
ters  more  than  chromaticity  values  per  se  in  deter¬ 
mining  whether  a  particular  purple  light  is  visually 
impressive. 

As  noted  above,  the  purple  light’s  increasing  satu¬ 
ration  is  easy  to  misinterpret  as  increasing  bright¬ 
ness.  Figure  6  shows  that  both  purple-light 
luminance  and  integrated  radiance  L  indeed  de¬ 
crease  steadily  with  ho-  Figure  5’s  chromaticity 
curve  reverses  once  again  at  ho  =  —7.19°,  presum¬ 
ably  marking  the  start  of  the  second  purple  light. 
The  chromaticity  curve’s  erratic  wiggles  here  are  in¬ 
structive:  The  weak  (and  thus  noisy)  signal  of  the 
ordinary  second  purple  light  is  barely  detectable  by 
instrument  and  certainly  not  by  naked  eye.  At  this 
ho,  skylight’s  mesopic-level  luminances  (~1  cd/m^) 
made  color  vision  difficult  and  chromaticity  moot: 
The  sky  simply  looked  dark  to  us. 

At  maximum  purity.  Fig.  5’s  Granada  late-spring 


448 


APPLIED  OPTICS  /  Vol.  42,  No.  3/20  January  2003 


Fig.  7.  Temporal  trends  in  purple-light  chromaticities  measured 
during  adjacent  evening  (5  September  1998)  and  morning  (6  Sep¬ 
tember  1998)  twilights  at  a  rural  site  near  Marion  Center,  Penn¬ 
sylvania.  Unlike  those  in  Fig.  5,  these  chromaticities  are  derived 
from  horizontal  irradiances. 


purple  light  is  perceptibly  redder  than  its  early- 
winter  counterpart.  Although  Fig.  5  does  not  gain¬ 
say  Neuberger’s  essentially  correct  statement  that 
the  purple  light  occurs  less  often  in  spring  than  in 
autumn,  12  it  does  show  that  spring  purple  lights  are 
no  less  vivid.  Equally  interesting  is  the  Granada 
chromaticity  cusp  at  =  —1.32°  that  marks  the 
purple  light’s  start.  Unlike  chromaticity  cusps  ob¬ 
served  in  angular  scans  along  sky  meridians^^  or 
across  the  rainbow, these  sudden  color  reversals  in 
the  purple  light  are  not  accompanied  by  any  abrupt 
change  in  radiance.  To  the  naked-eye  observer,  twi¬ 
light’s  slow  progress  means  that  such  temporal  sta¬ 
tionary  points  in  sky  color  go  unseen.  We  examine 
the  physical  reason  for  such  chromaticity  cusps  in 
Section  6  below. 

4.  Varietal  Purple  Lights 

The  ordinary  purple  light  often  dominates  the  color  of 
ambient  illumination  during  twilight,  and  so  its  col¬ 
orimetric  signature  not  surprisingly  appears  in  chro¬ 
maticities  derived  from  horizontal  irradiances  (Fig. 
7).  Except  for  their  displacement  toward  the  blue. 
Fig.  7’s  temporal  chromaticity  curves  from  Marion 
Center  look  quite  similar  to  their  radiance-based 
counterparts  in  Fig.  5.  However,  there  are  some  im¬ 
portant  differences. 

First,  the  9-5-98  curve  shows  no  hint  of  a  second 
purple  light.  In  other  words,  on  this  date  daylight 
shows  no  consistent  reddening  early  in  evening  nau¬ 
tical  twilight.  This  makes  sense  because  low  irradi- 
ance  levels  stopped  our  measurements  at  = 
-6.82°,  slightly  before  the  second  purple  light  begins 
on  most  evenings  (on  average,  at  Hq  ~  -7.4°;  see 
Table  2).  Furthermore,  the  single-scattering  geom¬ 
etry  of  nautical  twilight  means  that,  unlike  the  main 


Table  2.  Summary  Statistics  for  Purple-Light  Geometry, 
12-29-99  to  11-15-01“ 


Parameter 

Start 

K  (°) 

Peeik 

End 

Gamut  (^) 

Minimum 

-3.33 

-4.69 

-8.46 

0.03229 

Maximum 

-0.22 

-1.43 

-5.54 

0.08645 

Mean 

-1.411 

-3.893 

-7.369 

0.05249 

Median 

-1.16 

-4.15 

-7.4 

0.04918 

Standard  deviation 

0.927 

0.709 

0.559 

0.01247 

“These  statistics  were  calculated  over  35  evening  purple  lights 
with  a  colorimetrically  distinct  start,  peak  (i.e.,  maximum  purity), 
and  end;  spectral  radiances  were  measured  at  /r  =  20°  and  4>rei  = 
0°.  Normalized  chromaticity  gamuts  is  defined  in  Ref  22  and  is 
a  measure  of  the  fraction  of  the  CIE  1976  UCS  diagram  that  a 
chromaticity  curve  spans. 


purple  light,  the  second  purple  light  is  confined  to  the 
near-horizon  sky  and  thus  cannot  dominate  horizon¬ 
tal  irradiances.  Of  course,  during  ordinary  twilights 
the  near-horizon  sky  lacks  any  vivid  color  at  such 
small  Second,  our  analyses  show  that  during 

clear  civil  twilights  approximately  35%-40%  of  the 
horizontal  irradiance  at  the  surface  comes  from  sky¬ 
light  radiances  at  s  45°,  with  the  result  that  hor¬ 
izontal  irradiances’  chromaticities  are  strongly 
infiuenced  by  these  high-/i  radiances.  During 
evening  twilight,  the  sky  becomes  purplish  sooner  at 
higher  h  because,  at  any  given  moment,  a  fixed  sunlit 
altitude  (say,  the  tropopause)  has  a  local  Sun  eleva¬ 
tion  that  decreases  as  we  look  toward  higher  h.  So, 
as  twilight  progresses,  reddened  sunlight  illuminates 
the  lower  stratosphere  sooner  at  higher  h,  and  we 
first  see  the  purple  light  there.  Because  Fig.  7’s  pur¬ 
ple  lights  are  influenced  by  skylight  coming  from  high 
h,  they  start  before  the  purple  lights  of  Fig.  5.  Fi¬ 
nally,  even  though  Fig.  7’s  twilights  are  separated 
only  by  —10  h,  notice  that  their  chromaticity  curves 
differ  by  >1  JND  at  most  h^.  Such  dusk-to-dawn 
variability  must  be  considered  as  we  model  the  ordi¬ 
nary  purple  light. 

Figure  8  demonstrates  that  the  zenith  sky  itself 
can  be  purplish.  Here  we  measure  spectral  radi¬ 
ances  within  a  1°  FOV  centered  on  the  zenith,  a  di¬ 
rection  where  existing  accounts  disallow  the  purple 
light.  Even  at  its  most  vivid  (on  1-21-01),  this  zenith 
purple  light  has  a  differently  shaped  chromaticity 
curve  than  its  low-/i  kin.  As  expected,  the  zenith 
curves  are  mostly  bluer  than  those  shown  in  Fig.  5. 
However,  the  characteristic  purple-light  excursion 
across  the  Planckian  locus  remains  on  1-21-01,  and 
our  color  notes  judged  this  zenith  sky  to  be  “distinctly 
purplish”  when  =  -3°.  Yet  a  mere  two  days 
later,  on  1-23-01,  the  Owings  chromaticity  curve  is 
entirely  to  the  left  of  the  Planckian  locus  and  the 
zenith  sky  looked  only  blue. 

At  this  point,  readers  accustomed  to  colorimetric 
absolutism  may  be  reeling.  How  can  we  move  chro¬ 
maticity  curves  around  the  CIE  color  space,  all  the 
while  using  the  single  color  adjective  “purplish”  to 
describe  them?  In  fact,  our  colorimetric  relativism 
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Fig.  8.  Temporal  trends  in  evening  twilight  chromaticities  mea¬ 
sured  at  the  zenith  at  Owings,  Maryland,  on  2,  21,  and  23  January 
2001.  Although  the  21  January  twilight  has  a  distinct  purple 
light,  two  days  later  there  is  none. 


CIE  1976  u' 

Fig.  10.  Comparison  of  temporal  trends  in  evening  twilight  chro¬ 
maticities  at  Granada  on  7-8  June  2000  (h  =  20°,  =  0°).  Note 

that  the  vivid  purple  light  of  7  June  largely  disappears  by  the  next 
day. 


is  just  a  practical  example  of  color  constancy:  Even 
though  in  principle  our  measured  chromaticity 
curves  differ  perceptibly  from  day  to  day,  we  judge 
color  in  its  immediate  spatial  and  temporal  context. 
As  noted  above,  the  chromaticity  curve’s  shape  tends 
to  be  more  significant  visually  than  its  actual  chro¬ 
maticities.  Thus  different  twilight  chromaticities  on 
different  days  all  can  look  purplish. 

How  does  a  given  twilight’s  purple  light  change 
with  h?  Figure  9  answers  this  question  for  the 
evening  of  6-7-00  at  Granada.  To  create  Fig.  9  we 
simultaneously  measured  l°-FOV  radiances  at  h  = 
20°,  45°  with  our  two  PR-650  spectroradiometers  and 


CIE  1976  u' 

Fig.  9.  Temporal  trends  in  evening  purple-light  chromaticities  at 
Granada,  Spain,  on  7  June  2000  for  horizontal  daylight  irradiances 
and  for  skylight  radiances  at  h  =  20°,  46°  and  <|)rei  =  0°- 


horizontal  irradiances  with  a  third  spectroradiom- 
eter,  a  LI-COR  model  LI-1800. Figure  9’s  most 
striking  feature  is  the  close  similarity  between  the 
curves  from  horizontal  irradiances  and  from  radi¬ 
ances  at  h  =  45°  (the  former  curve  ends  at  /iq  = 
—3.36°  because  of  the  LI-1800  instrument’s  lower 
sensitivity).  This  similarity  reinforces  our  remark 
about  Fig.  7  that  chromaticities  at  high  h  infiuence 
the  color  of  horizontal  irradiances  (i.e.,  the  ambient 
color  of  twilight).  And,  although  Fig.  9’s  h  =  20° 
curve  is  much  redder  than  the  other  two,  all  three 
kinds  of  skylight  undergo  similar  changes  in  color 
during  twilight.  This  makes  good  physical  sense, 
because  all  skylight  derives  from  a  similar  sequence 
of  color  changes  undergone  by  the  setting  Sun  as  it 
illuminates  different  altitudes  in  the  atmosphere. 

Whatever  the  value  of  h,  marked  day-to-day  vari¬ 
ability  is  the  norm  for  the  ordinary  purple  light.  In 
Fig.  10,  two  twilight  chromaticity  curves  pair  the 
vivid  Granada  purple  light  of  6-7-00  with  a  nearly 
invisible  one  on  6-8-00  (the  6-7-00  curve  also  appears 
in  Figs.  5  and  9).  Even  though  both  days  were  cloud¬ 
less,  visibility  was  distinctly  better  and  atmospheric 
turbidity  lower  on  6-8-00.  Thus  one  visual  harbin¬ 
ger  of  a  pronounced  purple  light  may  be  a  daytime 
zenith  sky  whose  blue  looks  more  desaturated  than 
normal.  As  a  corollary,  any  twilight  chromaticity 
curve  with  a  cusp  near  sunset  or  sunrise  may  indicate 
the  presence  of  subvisual  cirrus  (see  Section  6  below). 

Figure  11  shows  just  how  variable  the  purple  light 
is,  even  when  it  is  visible.  In  Fig.  11  we  plot  with  X’s 
the  chromaticities  of  the  most  vivid  purples  in  35 
different  twilights  with  visible  purple  lights.  Collec¬ 
tively,  the  chromaticities  of  these  apparently  pur¬ 
plish  twilights  straddle  the  Planckian  locus  and  have 
a  root-mean-square  colorimetric  gamut  g  = 
0.01371,22  or  ~8.1  MacAdam  JNDs  at  this  u' ,  v'. 
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CIE  1976  u' 

Fig.  11.  Variability  of  35  purple-light  maxima  at  three  of  our 
ohserving  sites.  Complete  temporal  chromaticity  curves  are 
drawn  for  Owings,  Maryland,  on  2-3  May  2000,  two  adjacent  days 
with  the  largest  colorimetric  difference  in  purple-light  maxima 
measured  from  December  1999  to  November  2001. 


We  draw  chromaticity  curves  for  our  two  extreme 
purple-light  cases,  5-2-00  and  5-3-00. Although  in 
one  day  the  purple-light  maximum  underwent  the 
greatest  colorimetric  shift  that  we  measured  in 
nearly  two  years,  visually  this  distinct  color  shift  did 
not  seem  unusual.  Thus  Fig.  11  illustrates  a  funda¬ 
mental  truism  about  the  ordinary  purple  light: 
However  visually  impressive  a  given  purple  light  is, 
spectrally  based  colorimetry  may  reveal  a  quite  dif¬ 
ferent  twilight  from  the  one  that  you  think  you  saw. 

What  do  purple-light  spectra  themselves  look  like? 
Most  resemble  the  maximum-purity  example  shown 
in  Fig.  12,  measured  at  Owings  on  5-24-00  when  = 
—4.31°.  Colorimetrically  speaking,  this  is  a  proper 
purple  rather  than  some  exotic  metamer:  Com¬ 
pared  with  the  spectrum  of  sunlight  above  the  atmo¬ 
sphere,  the  purple-light  spectrum  has  less  energy  in 
the  yellows  and  yellow-greens  than  in  the  reds,  or¬ 
anges,  and  blues.  The  colorimetric  variability  seen 
in  Fig.  11  arises  largely  from  shifting  balances  be¬ 
tween  Fig.  12’s  two  broad  peaks  near  its  spectral 
extremes,  and  there  are  also  daily  fluctuations  in  the 
details  of  atmospheric  absorption  bands  at  wave¬ 
lengths  of  >660  nm. 

5.  The  Troposphere’s  Role  in  Causing  the  Purple  Light 

Figure  11  poses  a  problem  for  explanations  of  the 
ordinary  purple  light  that  rely  heavily  on  the  strato¬ 
sphere.  Although  we  cannot  rule  out  a  large  24-h 
increase  in  stratospheric  aerosol  from  2  to  3  May 
2000,  this  seems  an  improbable  first  solution.  Be¬ 
cause  Fig.  ll’s  chromaticity  extremes  occur  within 
24  h  and  because  there  is  no  record  of  extreme 
changes  in  the  stratosphere  during  that  time,  the 
troposphere  emerges  as  a  likelier  source  of  this  sud¬ 
den  shift  in  twilight  color. 


Fig.  12.  Visible-wavelength  spectral  radiances  for  a  typical 
purple-light  maximum  measured  at  Owings  on  24  May  2000  when 
ha  =  —4.31°  compared  with  those  of  sunlight  above  the  atmo¬ 
sphere  (normalized  to  have  the  same  sum  as  the  purple-light  spec¬ 
trum).  Of  the  35  purple-light  maxima  plotted  in  Fig.  11,  this 
purple-light  spectrum  has  the  chromaticity  that  is  closest  to  the 
mean  u' ,  v' . 


In  fact,  data  from  the  HALogen  Occultation  Exper¬ 
iment  (HALOE)  limb-sounding  satellite^^  show  that 
both  the  differential  profiles  and  the  total  column 
amounts  of  stratospheric  aerosols  at  our  site  can  be 
greater  when  there  is  no  purple  light.  Just  the  op¬ 
posite  would  hold  true  if  stratospheric  scattering 
were  the  sole  (or  even  the  chief)  factor  in  causing 
vivid  ordinary  purple  lights.  For  example.  Fig.  13 
shows  temporal  chromaticity  curves  for  three  differ¬ 
ent  twilights.  Although  the  USNA 11-12-98  twilight 


Fig.  13.  Temporal  trends  in  evening  twilight  chromaticities  at  the 
USNA  (12  November  1998)  and  at  Owings  (29  and  31  December 
1999).  S(i  is  the  total  column  amount  of  stratospheric  aerosol 
measured  above  these  sites  (see  Fig.  14  for  details). 
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Fig.  14.  HALOE  limb-sounding  satellite  measurements  of  the 
vertical  distribution  of  normalized  aerosol  surface  area  (i  (square 
micrometers  per  cubic  centimeter)  in  the  stratosphere,  interpo¬ 
lated  to  the  latitudes  and  longitudes  of  the  USNA  (on  12  November 
1998)  and  of  Owings  (on  29  and  31  December  1999).  Compare  the 
total  column  amounts  Sp  of  these  scatterers  with  twilight  colors  on 
the  same  dates  in  Fig.  13. 


has  no  purple  light,  the  Owings  12-29-99  purple  light 
(see  Fig.  5  for  a  closer  view)  is  vivid  but  clearly  less 
red  than  that  on  12-31-99. 

Now  compare  these  chromaticity  curves  with  Fig. 
14’s  HALOE  vertical  profiles  of  stratospheric  aerosols 
on  the  same  dates.  These  aerosols  are  primarily  sul¬ 
furic  acid  droplets,  and  their  median  radii  are  smaller 
than  0.1  |xm  at  all  altitudes  z.  Summing  these  drop¬ 
lets’  normalized  surface  areas  (3  (in  square  microme¬ 
ters  per  cubic  centimeter)  for  z  from  14.4  to  30.3  km 
(Ref.  30)  reveals  that  11-12-98  has  the  largest  column 
total  Sp  and  12-31-99  has  the  smallest  S(3.  This 
result  is  opposite  that  expected  if  stratospheric  aero¬ 
sol  concentration  alone  determined  the  vividness  (or 
even  the  existence)  of  the  ordinary  purple  light.  Nor 
does  the  vertical  distribution  |3(z)  solve  this  conun¬ 
drum:  When  the  12-31-99  purple  light  is  at  its  pur¬ 
est  (/iq - 4.46°),  all  sunlight  rays  that  pass  through 

the  troposphere  en  route  to  our  radiometer’s  LOS 
intersect  that  line  at  z  <  21  km.  These  rays  pass 
beneath  tropopause-grazing  ray  c  in  Fig.  4.  Yet  (3(z) 
below  21  km  usually  is  smaller  on  12-31-99  than  on 
Fig.  13’s  other  two  dates.  So  neither  the  vertical 
distribution  nor  the  column  totals  of  stratospheric 
aerosols  can  by  themselves  account  for  its  purple 
lights. 

Because  spectrally  selective  extinction  is  much 
greater  in  the  troposphere  than  in  the  stratosphere, 
the  troposphere  provides  the  reddest  illumination  of 
stratospheric  aerosols  along  our  LOS.  (We  support 
this  claim  quantitatively  in  Section  6  below).  If  the 
stratosphere  does  not  predominate  in  creating  ordi¬ 
nary  twilights’  colors,  then  the  troposphere  must. 
Figures  13  and  14  provide  indirect  evidence  of  the 


Fig.  15.  Photograph  of  a  purplish  twilight  crepuscular  ray  at 
Owings  on  1  January  2002  at  2222  UTC  (h^  =  —5.43°).  White 
lines  have  heen  drawn  on  the  sky  to  help  to  locate  this  low-contrast 
ray. 


troposphere’s  importance  in  influencing  twilight  col¬ 
ors,  including  the  purple  light.  Direct  visual  evi¬ 
dence  comes  from  Fig.  15’s  twilight  crepuscular  ray. 
Often  noted  by  purple-light  observers,^^  these  low- 
contrast  rays  can  span  the  entire  sky  as  alternating 
bands  of  pinkish-purple  and  blue:  i.e.,  the  purple 
light  flanked  by  its  absence.  Shadows  cast  on  the 
lower  stratosphere  eliminate  the  purple  light  within 
the  blue  bands,  and  mountains  and  tropospheric 
clouds  far  beyond  the  sunset  horizon  usually  cast 
these  shadows.  Thus  the  purple  light  disappears  if 
the  stratosphere  above  us  is  shaded  from  sunlight 
that  has  been  reddened  by  transmission  through  the 
troposphere. 

We  can  best  explain  Fig.  13’s  variable  purple  light  by 
examining  Fig.  16.  Each  of  these  geostationaiy  sat¬ 
ellite  images  was  taken  within  a  few  hours  of  our  cor¬ 
responding  twilight  measurements  in  Fig.  13.  In 
each  part  of  Fig.  16  we  have  drawn  a  red  arrow  from 
our  observing  site  toward  the  Sun’s  azimuth  when 
/iq  =  -4.5°,  i.e.,  toward  the  Sun  when  the  purple  light 
is  near  its  peak.  The  arrow’s  tip  marks  the  point 
where  the  Sun’s  unrefracted  rays  would  then  be  tan¬ 
gent  to  the  Earth’s  surface.  Notice  that  only  for  the 
two  December  1999  twilights  does  a  long  cloud-free 
path  extend  from  our  site  toward  sunset.  On  those 
days,  tropospherically  reddened  sunlight  can  illumi¬ 
nate  the  lower  stratosphere  along  our  LOS,  and  thus 
the  purple  light  forms.  But  on  11-12-98  [Fig.  16(c)]  a 
large  cirrus  shield  to  our  southwest  prevented  direct 
sunlight  from  being  transmitted  along  long  slant  paths 
through  the  troposphere  to  our  LOS.  Although  sun¬ 
light  that  has  been  slightly  yellowed  by  transmission 
through  the  stratosphere  illuminates  our  radiometer’s 
LOS,  no  tropospherically  reddened  sunlight  reaches  it. 
As  a  result,  no  purple  light  formed  on  11-12-98.  In 
our  research,  this  cause-and-effect  relationship  be¬ 
tween  the  purple  light  and  clouds  far  over  the  sunset- 
sunrise  horizon  is  quite  consistent:  Every  clear 
twilight  without  a  purple  light  has  some  cloud  cover  in 
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Fig.  16.  Visible-  and  infrared-wavelength  geostationary  satellite  images  of  the  United  States  (a),  (b)  on  the  afternoons  of  29  December 
1999  (visible  wavelengths)  and  31  December  1999  (visible),  respectively,  and  (c)  in  the  early  evening  of  12  November  1998  (infrared).  Red 
arrows  point  from  our  observing  sites  toward  the  Sun’s  azimuth  at  =  —4.5°.  Clouds  far  beyond  our  sunset  horizon  on  12  November 
eliminated  the  purple  light  then  (see  Fig.  13). 


the  Sun’s  direction.  This  cloud  cover  need  not  be  ex¬ 
tensive  and  may  be  hundreds  of  kilometers  distant, 
but  it  is  always  present. 

In  fairness,  Bullrich  states  explicitly^®  and  Deir- 
mendjian  implies^®  that  distant  clouds  matter  in  ex¬ 
plaining  the  purple  light.  But,  until  now,  the  purple 
light’s  disappearance  usually  has  been  regarded  as 
little  more  than  a  matter  of  local  cloud  climatology. 
In  fact,  a  thin  to  moderately  thick  overcast  does  not 
rule  out  the  purple  light’s  formation.  Although  one 
cannot  see  the  clear  sky  through  an  overcast,  our  mea¬ 
surements  show  that  the  purple  light’s  diffuse  illumi¬ 
nation  above  the  clouds  makes  some  overcast  twilights 
(briefly)  at  least  as  purple  as  their  clear-sky  counter¬ 
parts.  Yet  whether  the  day  is  clear  or  cloudy,  the 
purple  light  disappears  if  clouds  anywhere  along  the 
Sun’s  azimuth  prevent  tropospherically  reddened  sun¬ 
light  from  illuminating  the  stratosphere  above  us. 

6.  Modeling  the  Purple  Light 

To  dissect  the  purple  light’s  complexities,  we  turn  to 
LOWTRAN  7,  a  radiative-transfer  model  that  incorpo¬ 
rates  considerable  meteorological  detail.®^'®®  To 
make  realistic  comparisons  with  Fig.  13’s  measure¬ 


ments  of  ordinary  purple  lights,  our  lowtran  simula¬ 
tions  in  Fig.  17  include  the  following  parameters:  a 
default  midlatitude  winter  atmospheric  profile  of 
temperature,  pressure,  humidity,  and  gas  mixing  ra¬ 
tios;  tropospheric  aerosols  typical  of  rural  sites;  back¬ 
ground  stratospheric  dust  and  other  aerosols  (i.e., 
ordinary  rather  than  volcanic  twilights);  no  clouds  or 
rain;  multiple  scattering;  a  surface  albedo  of  0.2;  and 
Mie  aerosol  phase  functions.  Consistent  with  our 
measurements  of  radiance  spectra,  the  model  is  set  to 
calculate  skylight  spectral  radiances  at  h  =  20°, 

4'rel  =  0°- 

Keeping  these  parameters  constant  in  Fig.  17,  we 
vary  Hq  and  horizontal  meteorological  range  As 
tropospheric  turbidity  and  mean  extinction  coeffi¬ 
cient  increase,  V  decreases  and  sunsets  tend  to  be¬ 
come  redder.  A  corollary  result  is  that,  as  V 
decreases,  irradiances  illuminating  stratospheric 
scatterers  along  our  LOS  during  twilight  are  redder 
at  a  given  Iiq.  Our  reference  case  in  Fig.  17  sets  V  = 
23  km,  a  value  that  is  typical  of  clear  rural  atmo¬ 
spheres.  All  the  features  seen  in  Fig.  13  are  repro¬ 
duced  in  this  case:  (1)  pre-sunset  reddening  is 
followed  by  rapid  bluing  of  skylight  before  sunset;  (2) 
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Fig.  17.  Radiative  transfer  model  lowtran  7’s  calculations  of 
temporal  chromaticity  trends  for  ordinary  purple  light  as  functions 
of  /iq  and  meteorological  range  V  at  h  =  20°  and  =  0°.  For 
comparison,  we  also  show  a  chromaticity  curve  calculated  for  a 
volccmic  twilight  with  high  levels  of  fresh  volcanic  scatterers  in  the 
stratosphere.  Note  the  similarities  between  the  F  =  23  km  curve 
(our  LOWTRAN  reference  case)  and  the  purple-light  chromaticities 
measured  in  Figs.  5  and  13. 


the  purple  light  starts  shortly  after  sunset;  (3) 
reaches  its  maximum  at  /iq  =  “4°,  and  (4)  is  followed 
by  a  second  purple  light  that  begins  near  civil  twi¬ 
light’s  end.  Decreasing  V  makes  the  purple  light 
redder  and  more  vivid  and  expands  its  gamut  (see  the 
y  =  5,  10  km  curves),  whereas  increasing  V  does  the 
opposite  (y  =  125  km  curve).  Although  these  re¬ 
sults  make  good  physical  sense,  a  less  intuitively  ob¬ 
vious  result  is  that  combining  background 
stratospheric  aerosols  with  a  turbid  troposphere  (y  = 
5  km)  yields  twilight  purples  that  match  the  most 
vivid  ones  seen  in  a  post-volcanic  twilight. 

We  have  claimed  that  tropospheric  extinction  pro¬ 
duces  most  of  the  red  component  seen  in  ordinary 
purple  lights  (see  Fig.  12).  Figure  18’s  lowtran  sim¬ 
ulation  provides  strong  colorimetric  support  for  this 
claim.  In  Fig.  18  we  have  redrawn  the  Planckian 
locus  and  Fig.  17’s  reference  case  at  a  much  smaller 
scale  to  show  the  differential  and  integral  colors  pre¬ 
dicted  as  we  move  down  along  Fig.  4’s  LOS  for  = 
—4°  (lowtran’s  purple-light  maximum).  The  stron¬ 
gest  source  of  skylight  radiances  in  our  LOS  is  singly 
scattered  direct  sunlight.  In  turn,  this  sunlight  has 
been  transmitted  through  various  optical  thicknesses 
of  the  atmosphere  to  reach  our  LOS.  We  calculate 
skylight’s  integral  colors  by  summing  net  skylight 
radiances  along  this  slant  path  from  the  top  of  the 
LOWTRAN  atmosphere  (100  km)  down  to  altitude  z  km. 
A  similar  calculation  yields  differential  colors  from  z 
to  z  +  1  km.  When  z  =  0  km,  the  model  by  definition 
reproduces  the  chromaticity  that  a  surface-based  ob¬ 
server  sees,  and  so  the  path-integrated  chromaticity 
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Fig.  18.  Differential  and  path-integrated  skylight  colors  pre¬ 
dicted  by  LOWTRAN  7  for  Fig.  17’s  reference  case  at  =  —4°  and 
along  the  line  of  sight  h  =  20°  and  <|)rei  =  0°-  In  this  figure  only, 
the  pluses  at  which  the  arrows  point  denote  particular  tangent-ray 
altitudes.  For  comparison,  we  also  show  the  entire  surface-based 
twilight  curve  for  the  reference  case  shown  in  Fig.  17. 


curve  intersects  the  reference  case’s  curve  (/iq  is  not 
fixed  for  the  latter). 

What  does  Fig.  18  reveal?  As  we  move  toward  the 
surface  along  our  LOS,  we  use  as  landmarks  the 
tangent  (i.e.,  lowest)  heights  of  sunlight  rays  that 
intersect  it;  the  intersections  themselves  occur  at 
higher  altitudes.  Beginning  with  the  ray  whose  re¬ 
fracted  tangent  height  z  =  65  km,  we  find  initially 
that  differential  skylight  colors  get  bluer  than  is  pos¬ 
sible  by  Rayleigh  single  scattering  (whose  chromatic¬ 
ity  corresponds  to  infinite  color  temperature).  This 
enhanced  bluing  apparently  is  caused  by  strato¬ 
spheric  ozone  absorption  in  the  spectral  Chappuis 
bands  centered  near  600  nm.®®  As  tangent  height  z 
decreases  still  further,  differential  contributions  to 
path-integrated  sky  color  reverse  near  u'  =  0.16;  they 
next  grow  greener,  and  then  yellower.  Ray  c  in  Fig. 
18  has  its  tangent  height  at  the  tropopause  (z  =  10 
km  in  lowtran;  Fig.  4),  and  it  is  the  last  stratosphere- 
only  contribution  to  the  path-integrated  skylight 
color  that  we  see  at  the  surface.  Rays  below  c  are 
transmitted  though  ever-increasing  optical  thick¬ 
nesses  of  the  troposphere. 

Although  the  differential  contribution  of  ray  c  is 
yellowish-green,  path-integrated  color  at  the  tropo¬ 
pause  along  our  surface-based  LOS  looks  distinctly 
blue.  Only  as  we  add  scattering  from  sunlight  rays 
that  have  traversed  the  troposphere  do  we  get  red 
contributions  to  surface  skylight.  Of  course,  be¬ 
cause  tropospheric  extinction  increases  rapidly  with 
decreasing  altitude,  tropospheric  rays  make  ever- 
redder,  but  also  ever-weaker,  additions  to  the  surface 
signal.  Yet,  even  after  scattering  from  Fig.  18’s  red¬ 
dish  surface-grazing  ray  has  been  added,  the  integral 
skylight  color  has  covered  only  approximately  half  of 
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Fig.  19.  LOWTRAN  7  calculations  of  the  dependence  of  the 
stratospheric  and  tropospheric  components  of  skylight  radiances 
received  at  the  Earth’s  surface.  The  model  atmosphere  here  is 
Fig.  17’s  reference  case,  and  the  observer’s  line  of  sight  is  h  =  20° 
and  4>rei  =  0°- 


the  colorimetric  distance  from  “K  to  its  final  purplish 
color  at  our  eye.  By  definition,  the  remaining  radi¬ 
ance  that  reaches  us  must  come  from  multiple  scat¬ 
tering  within  the  Earth’s  shadow,  and  this  last 
differential  contribution  is  purplish  to  blue  (Fig.  18, 
lower  left).  The  stratosphere  is  far  from  irrelevant 
here,  because  it  provides  some  of  the  bluish 
molecular-scattering  contribution  seen  at  Fig.  12’s 
short  wavelengths.  Yet  lowtran  indicates  that  scat¬ 
tering  and  extinction  within  the  troposphere  (includ¬ 
ing  multiple  scattering)  are  by  far  the  most  important 
processes  in  creating  the  ordinary  purple  light’s  red 
components. 

Finally,  what  can  we  tell  from  lowtran  about  the 
chromaticity  cusps  seen  in  Figs.  5  and  7?  Part  of  the 
answer  comes  in  Fig.  19,  where  we  plot  as  functions 
of  ho  two  parts  of  the  spectrally  integrated  lowtran 
radiances  received  at  the  surface.  We  call  the  part 
of  path-integrated  L  that  originates  solely  in  the  tro¬ 
posphere  <  10  km).  The  other  part  is  > 
10  km),  and  it  originates  in  the  stratosphere;  its  Fig. 
19  values  refiect  extinction  losses  after  transmission 
through  the  troposphere  to  the  surface.  Thus  Fig. 
19  tells  a  simple  story:  At  ho  ~  -0.65°,  the  strato¬ 
sphere’s  contributions  to  surface  radiances  first  ex¬ 
ceed  those  from  the  troposphere.  At  still  lower  ho, 
stratospheric  radiances  begin  to  outweigh  their  tro¬ 
pospheric  kin. 

Although  this  result  may  seem  to  contradict  Fig. 
18’s  arguments,  in  fact  it  does  not:  For  Fig.  18  it 
made  sense  to  discuss  tangent  heights  for  sunlight 
rays  that  intersect  our  LOS.  At  ho  =  -4°,  these 
tangent  heights  are  much  lower  than  the  subsequent 
LOS  intersection  heights,  which  is  the  altitude  given 
in  Fig.  19.  For  example,  when  h  =  20°  and  ho  = 
-4°,  the  point  whose  z  =  10  km  on  our  LOS  is  actu¬ 
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Fig.  20.  LOWTRAN  7  calculations  of  the  effects  of  subvisual  cirrus 
on  purple-light  chromaticity  trends,  including  the  formation  of  a 
chromaticity  cusp  near  sunset.  The  uniform  cirrus  layer’s  base  is 
at  10  km,  and  its  normal  optical  thickness,  t,  is  0.02. 


ally  several  kilometers  within  the  Earth’s  geometric 
shadow  (i.e.,  this  point’s  “tangent  height”  is  negative). 

Figure  20  completes  the  important  puzzle  of  the 
chromaticity  cusps.  There  we  repeat  our  lowtran 
reference  case’s  chromaticity  curve  for  the  ordinary 
purple  light  and  label  it  “background  aerosol  only.” 
Next  we  add  a  500-m  thick  veil  of  subvisual  cirrus  at 
z  =  10  km  and  make  this  cloud’s  normal  optical  thick¬ 
ness  T  =  0.02,  a  realistically  small  value.  The  effect 
is  striking:  We  redden  the  start  of  evening  twilight 
and  considerably  sharpen  the  chromaticity  corner 
that  marks  the  purple  light’s  onset.  If  we  double  t  to 
0.04  (not  shown),  this  corner  changes  to  a  dagger- 
sharp  cusp,  a  likely  colorimetric  signature  for  subvi¬ 
sual  cirrus.  Why  does  this  sudden  transition  occur? 
We  infer  from  Fig.  19  that  stratospheric  and  tropo¬ 
spheric  contributions  to  surface  skylight  are  approx¬ 
imately  equal  at  sunset.  Now  we  have  added  a  thin 
layer  of  cirrus  at  the  tropopause.  Although  this  lay¬ 
er’s  Sun  elevation  is  slightly  higher  than  that  at  the 
surface,  its  forward  scattering  of  reddened,  grazing- 
incidence  sunlight  nonetheless  increases  sharply 
close  to  the  time  of  surface  sunset.  We  measure  the 
net  result  as  an  unusually  sharp  reversal  of  skylight’s 
pre-sunset  color  trend — and  as  an  abrupt  start  of  the 
purple  light.  Thus,  although  the  ordinary  purple 
light  may  never  be  a  mainstay  remote-sensing  tool,  it 
continues  to  present  us  with  an  unexpected  combina¬ 
tion  of  celestial  beauty  and  optical  insight. 

7.  Conclusions 

What  simple  physical  explanation  can  we  offer  for  the 
sequences  of  ordinary  twilight  colors  shown  in  Figs.  5 
and  7?  As  the  Sun  sinks  below  ho  =  5°,  rapid  in¬ 
creases  in  slant-path  optical  thickness  for  direct  sun¬ 
light  at  first  redden  the  irradiances  incident  on 
tropospheric  scatterers  (Fig.  21).  Although  sunlight 
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Fig.  21.  Temporal  trends  in  the  colors  of  twilight  irradiances 
illuminating  a  vertical  plane  that  faces  the  Sun,  measured  at  the 
USNA  on  21  September  and  12  October  2000.  When  <  l°-2°, 
the  irradiance  from  skylight  begins  to  exceed  that  from  direct 
sunlight,  thus  causing  scattering  by  tropospheric  scatterers  to  be¬ 
come  increasingly  bluish  then. 


radiances  grow  ever  redder,  their  rapidly  decreasing 
magnitude  means  that,  when  surface  =  l°-2°,  the 
incident  irradiance  on  these  scatterers  starts  to  be¬ 
come  bluer  (i.e.,  skylight  now  begins  to  dominate  this 
irradiance).  Soon  after  sunset,  the  still-illuminated 
stratosphere  begins  to  predominate  over  a  tropo¬ 
sphere  that  is  passing  rapidly  into  shadow  (Fig.  19), 
and  reddened  sunlight  sent  to  the  surface  by  strato¬ 
spheric  scatterers  results  in  the  reddening  of  skylight 
that  marks  the  purple  light’s  onset. 

Although  skylight  radiances  continually  decrease 
during  evening  twilight,  a  combination  of  (1)  tropo- 
spherically  reddened  sunlight  that  is  singly  scattered 
in  the  stratosphere  and  (2)  purplish  skylight  that  is 
multiply  scattered  in  the  troposphere  at  first  in¬ 
creases  this  skylight’s  red  component.  The  detailed 
vertical  distribution  of  stratospheric  scatterers  (pri¬ 
marily  sulfuric  acid  droplets)  appears  to  be  less  im¬ 
portant  in  causing  this  reddening  than  the  simple 
fact  that  the  scatterers  are  present.  For  Hq  <  -4°, 
process  (1)  above  fades  into  insignificance  as  the 
lower,  denser  parts  of  the  stratosphere  pass  into 
shadow.  Now  ozone  absorption  and  molecular  sin¬ 
gle  scattering  in  the  stratosphere  combine  with  blu¬ 
ish  multiple  scattering  in  the  troposphere  to  make 
skylight  increasingly  blue.  This  process  continues 
until  tropospherically  reddened  sunlight  illuminates 
a  stratum  of  scatterers  at  even  higher  altitudes.  For 
nonvolcanic  twilights,  however,  this  second  purple 
light  has  such  low  luminance  that  we  cannot  see 
either  its  color  or  its  enhanced  brightness. 

Thus  the  ordinary  purple  light  is  both  simpler  and 
more  complex  than  previously  thought:  Molecular 
scattering  of  sunlight  at  high  altitudes,  combined 
with  reddened  illumination  of  the  stratosphere  and 


the  troposphere,  produces  the  subtle  purples  of  non¬ 
volcanic  twilights.  LOWTRAN  7  shows  that  although 
stratospheric  scatterers  are  needed  to  send  reddened 
sunlight  to  the  surface,  at  most  they  yellow  (rather 
than  redden)  sunlight  transmitted  through  the 
stratosphere  alone.  The  resulting  purples’  details 
depend  largely  on  daily,  or  even  hourly,  changes  in 
the  troposphere’s  optical  properties  beyond  the  sun¬ 
set  horizon. 
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